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TWO ALTERNATIVE WAYS FOR SOLVING THE COORDINATION PROBLEM 
IN MULTILEVEL OPTIMIZATION 

Jaroslaw Sobicszczanski Sobicski 

August 1991 

Summary 

The paper describes two new techniques for formulating the coupling between levels in multilevel opti- 
mization by linear decomposition, proposed as improvements over the original formulation, now several years 
old, that relied on explicit equality constraints which were shown by application experience as occasionally 
causing numerical difficulties. The two new techniques represent the coupling without using explicit equality 
constraints, thus avoiding the above difficulties and also reducing computational cost of the procedure. The 
old and the new formulations arc presented in detail, illustrated by an example of a structural optimization. 
A generic version of the improved algorithm is also developed for applications to multidisciplinary system not 
limited to structures. 

Notation 

A vector of cross* sectional areas*, Ay 

Cy cumulative constraint of ith beam 

DISi vector defined by cq. (22)* 

g k vector of constraints* for a beam, c.g., stress limits, and local buckling, k = 1... NGB ) 

partitioned in subsets of lengths NGBy ) each subset corresponding to ith beam 


Gjs vector of constraints* for the assembled structure, c.g., displacement limits, k = 1 . . . NGA 

1 vector of cross-sectional moments of inertia*, Iy 

Li length of ith beam 

NE number of beams in a framework 

Ny vector of the end forces* for ith beam 

NPy length of vector Py 

NSS number of subsystems 

NX length of vector X 

Py vector of parameters* in optimization of ith beam, comprising elements P q , q = 1 . . . NPy 

qth element of vector Py 
SA system analysis 

SI input vector of length N SI into S' A* 

SO output vector of length NSO from SA* 

SSAy ith subsystem black box analysis 

{SSC} vector of geometrical* variables determining the structure overall shape 
SSIy input vector of length NSSIy into SSAi* 

SSOy output vector of length NSSOy from SSAy * 

TOL tolerance parameter set by user 



w 

X 


weight, equivalent to volume in a homogeneous structure 

vector of design variables* at the system level (assembled structure) in optimization bv 
decomposition J 


Yi vector of design variables* in tth beam optimization problem 

vector of design variables* in optimization without decomposition 
{ } brackets identifying vectors arc omitted where possible without causing ambiguity. 

dcfincTwhcrr^^ 2 ’ C S '’ XU ' XL ' d “° lC “ PPCr “ d '° WCr b ° u,lds ltcsc va ™blc3; other symbols ore 

Introduction 

fnn ^ gC ; scalc optimization problems benefit from decomposition into a set of smaller, more manageable 
w^rthTsv^ bl SU Pr °^ CmS ‘ i“ a u hicrarchic dccom Position method, the subproblems form a pyramid 
arc couoled th T °V°? SUbsyStCm problcms in thc horizontal layers below. These subproblems 
dceotnnrwV thr ° llgh thc solut,on of a coordination problem. A particular procedure for optimisation by 
dwomposition introduced in reference 1 solved thc coordination problem by enforcing a set of equality 

referent Zl ,CV f S “ d * llsi "S a " °P^um sensitivity analysis formulated in 

fnrmnU^ r r i T pr f ° ccdurc > rcfcrrcd to “ Optimization by Linear Decomposition (OLD), was 
formulated for two-levels in reference 1 and was demonstrated using a framework as a test case in reference 4 
representing a class of skeletal, redundant structures. It was subsequently generalized to an arbitrary number 
cvcls in reference 5. The OLD is a generic method applicable to any system that is amenable to a hierarchic 
decomposition, c.g., multidisciplinary applications reported in references 6 and 7. 

Practical experience with thc procedure, and its examination in reference 8 point to the enforcement of 

Snvcree'Th m i S “ * S ° UrCC ° f “ Umcrical difficu,tics that occasionally make this procedure slow to 
converge. This observation is consistent with thc opinion generally held among the developers and users 

o optimization methods that addition of equality constraints to an optimization problem tends to make thc 
solution numerically more difficult. 

infr^° tlV ? Cd b r tHC ab ° VC ’ bytW ° nCW tcchni( l ucs satisfying equality constraints that were recently 
duccd in references 9 and 10, this paper defines two alternative modifications to thc OLD procedure Either 

riplSr rCm ° VCS CXP C<1Uality COnStraint8 from thc P^ecdure while, still, achieves their enforcement 

The two alternative techniques will be introduced by generalizing from an example of the framework 

— h iT a T° “ rcfcrcnccs 1 and 4 ' To this cnd > thc framework analysis will be discussed first, 
followed by the framework optimization problem formulated in a standard manner without decomposition 

Next, an abridged description of thc OLD will be given, limited to a twolcvcl, structural optimization of the 
introduced ^ ^ Wlth ^ “ & rcfcrcncc> ^ ncw alternative modifications to thc OLD algorithm will be 

Original Reference Procedure 
Analysis 

The framework is shown in figure 1. We limit thc framework analysis to a twodimensional case by allowing 
only in-plane displacements under thc action of static loads. Thc static analysis of this structure may be 
ormulatcd in two levels by using either a substructuring approach where each of the three beams in the 
framework is regarded as a separate substructure, or a finite element method in which each beam is a single 
finite element. Choosing the finite element formulation, the analysis of thc assembled framework requires an 
input and generates an output as defined in table 1. 
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ni r 

The derivatives such as the displacement U with respect to the cross-sectional area A j, ^7, exist because 
l) = /( A , /, A'.STy), /Vj = /(/l, I,SSC), and L( = /(55(7). These derivatives arc obtained by finite differencing 
or by a quasi-analytical sensitivity analysis embedded in the framework analysis. 

Selected data from the above output arc entered as input into a local strength analysis of the ith beam. 
The level of refinement of that analysis is immaterial for the purposes of this discussion, c.g., it might be a 
finite element analysis whereby the beam is divided into number of finite elements, or an elementary strength 
of materials analysis. In either ease, the input and output arc as defined in table 2. 

The above approach suggests decomposition of the framework analysis into the assembled framework 
analysis (system level analysis) and the beam analyses (subsystem level analyses). This decomposition forms 
a hierarchy in which the former is a parent and the latter arc the daughters. The parent-daughter relationship 
is hierarchic because the information flows from the parent to each daughter and no information is directly 
transmitted from one daughter to another. Since the daughter analyses are mutually independent, they may 
be executed concurrently. 

To make the system and subsystem level analyses consistent, one has to acknowledge that A{ and are 
functions of the beam cross-sectional dimensions {V r }j. Therefore, when prescribing the values for -A*,/*, and 
{Y)i, we must satisfy the consistency equations 

Ai=fA({y)ih h = fdiY)i ), i = l...NE ( 1 ) 

Keeping these equations satisfied in a multi-level optimization is a part of the so-called coordination problem, 
and the analysis by decomposition will be exploited to establish a corresponding optimization-by-decomposition 
scheme that will be described later. 

Optimization Without Decomposition 

The optimization problem to be solved is 

“find Z such that F(Z) is at minimum subject to constraints” (2) 

Gj(Z)< 0, j = l...NGA; (2a) 

g k (Z)< 0, A: = 1 . . . NGB (2b) 

ZL < Z < ZU (2c) 

In the static problem at hand, the design variables may include the elements of the {SSC} and some, or 
all, of the cross-sectional dimensions of the beams, {y}*. The constraints G * comprise the assembled structure 
displacements and clastic stability, and at the beam cross-section level the constraints entail the allowable 
stress, beam column buckling, and local buckling. 

Optimization by Linear Decomposition (OLD) 

The above problem is decomposed into a single problem solved for the assembled structure (system level 
problem) and NE separate problems solved for each beam (subsystem level problems). It is convenient to 
describe the subsystem level problem first. The equation underscoring indicates which parts of the formulation 
will be changed by the two new techniques formulated in this paper. 

Subsystem level . 

At the subsystem level, an optimization problem for the ith beam is independent of the other problems at 
that level, hence the beam optimizations may be executed concurrently. The ith beam problem is, taking into 
account cq. (1): 


“find {y},* such that Ci({Y}i } {P}i) is at minimum subject to constraints” (3) 
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Ai -f A ({Y)i) = 0 , ( 3a ) 

I J -fr{Y}i = 0, (3b) 

{YL)i < {YU < {YU} i; ( 3c) 

where {YL}i and {YU)i arc the lower and upper side-constraints on {y}<. The {/>}< is a vector of parameters 
comprising A { J„ L u and {N}i that arc output from the analysis at the system level as defined in table 1 
and arc passed as input into the beam analysis as shown in table 2. These parameters stay constant in the 
process of the beam optimization. Hence, for tth beam 

{Pq}i = {A„/ f , L it {AT}*}, q = l...NPi (4) 

The scaiar Q is a cumulative constraint representing the degree of satisfaction, or violation, of the subset 
9k(i /t'/> * - 1 .. . NGBi , pertaining to the tth beam, of the entire set of g k that appears in cq. (2b). The 

cumulative constraint may be evaluated as in references 4 and 5 by means of the Kreissclmeier-Steinhauser 
(KS) function (ref. 11): 


C d9k) = KS(g k ) = (1/p) In exp (pg k )j , k = 1 . . . NGB { 


( 5 ) 


where p is a user-controlled coefficient that governs the distance between the KS and max(^). An alternative 
K 5-formulation that avoids generation of large values of the exponential function is 


c i{9k) - maxfe) + 1/p) In ^J^exp (p (g k - max(g k )))^j , k = 1 . . . NGB { 
The KS function is differentiable and has the property that 


(5a) 


max(p fc ) < KS(g k ) < max^*) + In NGB { /p, k = 1 . . . NGBi (6) 

approximating the nondiffcrcntiablc max(« 7 *) with an error dependent on p (the larger p is, the smaller the 
error is. However, optimization may be more difficult numerically). 

The optimization in cq. (3) alters {Y}i to enforce the consistency cq. (1) by means of the equality constraints 
in cqs. (3a) and (3b), and to achieve a minimum of C, equivalent to a minimum ofmaxfa*). This optimization 
produces an optimal solution comprising of C imin and { Y } iopt . It is followed by an optimum sensitivity analysis 

(ref. 2) to obtain the derivatives of Cj min and {Y } iopt with respect to {Pq)i, denoted and 

r uPqi 

System level. 

r /oo^ 87810111 lcvcl P roblcm > thc design variable vector X contains Ai,I { , i = 1 . . . NE, and the elements 
or (ojC } 

{A} = {A i ,/i,{55C}}, i = l...NE-, ( 7 ) 

At this level, the assembled structural analysis whose input is defined in table 1 is carried out. Thc objective 

function in this ease is thc structural weight which for a homogeneous material may be replaced with the 
material volume 

NE 

F(X) = W(X) = '£A i L i (8) 
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The problem formulation in 


“find X such that F(X) is at minimum subject to constraints” (9) 

G k (X)<0, k = l...NGA; (9a) 

Cj{ X) < TOL, i = l...NE ; (9b) 

{YI'}i< {Y)i<{YU}i, i = l...NE ; (9c) 

XL < X < XU ; (9d) 

where TOL is a suitable tolerance parameter. 


The constraints G^X) pertain to the assembled structure behavior, c.g., displacement limits and overall 
clastic stability. The C t (X) is approximated by extrapolation 

c,m = c imin +-£T. - x j0 ), 

i=l...NE, j = l...NX, q=l...NPi (10) 

where P qi is defined by cq. (4). In the above, C tmin and - g) P ^ arc transmitted from the subsystem level 
optimization and optimum sensitivity analysis, and the derivative product terms constitute a chain rule 
differentiation necessary since some of the parameters P q depend on Xj as noted in discussion of table 1. 

Since Xj is an element of {P}*, we will have 3-^7 = 1 for the coincidences (P 9 )j = Xj that occur in the 

summation. The values of ( Xj) 0 arc those for which the assembled framework analysis was carried out prior 
to the current optimization. 

When C x is expressed by a KS function as in cq. (5) or (5a), the extrapolation error in eq. (10) may be 
significantly reduced, or shown in references 4 and 5, by extrapolating each constraint function g k that enter 

Q 


9k(X) = gko + EE «-**>■ 

J = 1 - - • NX, q = l...NPi (10a) 

and, then, computing in extrapolated C x by means of either cqs. (5) or (5a). This technique removes that part 
of the extrapolation error that in cq. (10) would be caused by the curvature of the logarithm and exponential 
functions embedded in the KS function. 

An extrapolation similar to cq. (10) is used to approximate {T}j in cq. (9c): 

xm = y jopt + 

i = l...NE, j = l...NX t q = l...NPi (11) 

XL<X < XU 

Finally, XL and XU, arc, respectively, the lower and upper limits on X. These side constraints include the 
move limits guarding against excessive extrapolation errors in cqs. (10) and (11). 
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Overall procedure. 

The overall OLD procedure is 

1. Initialize the overall shape data in {SSC}, and the beam cross-sectional dimensions (y} i; 

2. Enter {T}* into cqs. (3a) and (3b) to initialize A and /; 

3. Execute the assembled structure analysis (table 1); 

4. For each beam execute a subsystem level optimization per cq. (3), that refers to the beam analysis per 
table 2, and carry out the optimum sensitivity analysis with respect to parameters defined in eq. (4); 

5. Execute optimization at the system level per cq. (9); 

6. If termination criteria set by user arc not satisfied, reset X and Y to the new values and repeat from 

step 3, otherwise exit. y 

To conclude the description of the OLD procedure at two levels, one should point out that the individual 
optimizations at both levels arc coupled by means of cqs. (3a), (3b), (9b), (9c), (10), and (11). These equations 
represent the coordination problem which is solved by virtue of converging the overall procedure. 

Shortcomings of the OLD Procedure in Need for Improvement. 

In the subsystem level formulation in eq. (3), there is a possibility of a conflict between the equality 

Tfe^ble s : t CqS ' 1 ( I’ and th f SldC Conatraints in «»• ( 3c )- Specifically, it may not be possible^ find 
collet the 1 I 0 " 1 ° .V T! 0 " 1 u Pr ° bICm WhUc Sati3fying both scts of constraints. To alleviate that 
eaf ?cl Ld tflW 1 T abov 7 ncludc8 approximate representation of the side constraints on {F}„ 

the A Tr 1 1 ’ °^ CCP hC Sy ? CVC °P tlndza * ion from imposing on the ith beam such combinations of 
the Ai and 17, values that cannot be attained with the physically realizable {y}<. The above potential conflict 
is one disadvantage of using the equality constraints in cqs. (3a) and (3b). 

As . mcntloncd m Introduction, the other disadvantage is the increased difficulty of solving the 
optimization problem of cq. (3) brought about by the presence of the equality constraints, as pointed out in 

constraints ^ ^ modificatio °s arc introduced next, primarily, to remove these equality 

Proposed Modifications 

Wh ° SC introduction is this P a P cr Purpose arc defined herein. The common 
feature of both modifications is removal of the equality constraints in cqs. (3a), (3b), (9b), and (9c) and 

an indirect fulfillment of these constraints by reformulating the subsystem optimization problem. The two 
modifications differ in the details of that reformulation. 

Modification 1 

The first modification is based on a technique for locating simultaneous roots of a set of functions using the 
KS function as described m reference 9. Specifically, if a set of NF functions FAY) = 0 for Y = Y 0 then a 
KS function comprising the positive and negative F^Y) has a minimum at Y a Formally 


If Fi(Y) = 0, for y = Y 0 , then KS(Fi(Y) - Fj(y)) 
is at minimum for Y = Y 0 \ i = 1 . . . NF 


( 12 ) 


The above property of the KS function may be used to satisfy the equality constraints in eqs. (3a) and (3b). 
ror brevity, we define nondimcnsional functions v ' 
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Fm = (Ai - f A ({Y}i))/Ai, i = l... NE 
Fn = (/< - fi({Y}i))/Ii, t = l... NE 


( 13 ) 



and construct a composite function 


c i ({V}<, {P}i) = KS(F Ai - F Ai ,F H - F,i) 


= (“) ln ^ CXP ( pFA d + cxp (~ pF/ *») + CX P (P F a) + CX P i-pFli)) 


(14) 


By virtue of cq. (12), this function is a minimum at {Y 0 }i where F/j and F/, also vanish. We will seek 
{Yo)i as {Y} iopt in an optimization that entails Ci({Y}j, {P}i) defined by cq. (14) as the objective function. 
This approach was shown to be effective in reference 12, although the details of the function formulation and 
of the overall procedure defined in that reference were different. 

The subsystem level optimization of cq. (3) for the ith beam will now change to the following one in which 
C{ is defined by cq. (14): 

“find {Y}* such that ^({Yjj, {F}*) is at minimum subject to constraints” (15) 

9k({y}i)< 0, k = l...NGBi (15a) 

[YL} X < {Y} { < {YU)i (15b) 

The subsystem level optimization satisfies the local constraints, cqs. (15a) and (15b), and comes as close to 
Ffa — 0 and Fj x = 0, that is fj = /*, as possible. However, in contrast to cqs. (3a) and (3b), it is not 

required to nullify these quantities completely, if that is not yet possible in the process of iterating between 
the levels in the overall procedure. Hence, the potential conflict among the constraints cqs. (3a), (3b), and 
(3c) described in the preceding section has been removed. Consequently, the constraints in eq. (9c) that were 
needed to alleviate that conflict are now deleted from the system level optimization. As a result of the deletion 

of cq. (9c) from the system level optimization, the derivatives become unnecessary. Consequently, 

the optimum sensitivity algorithm from reference 2 may be replaced with the computationally less expensive 
algorithm from reference 3 for a significant reduction of the computational cost of the entire procedure. 

The system level optimization change is limited to deletion of cq. (9c), and to redefinition of C{ from the 
one given in cq. (5) to that set by eq. (14). The new definition of C x must also be used in eq. (10). 

Even though the equality constraint such as those in eq. (3) do not appear directly in the above subsystem 
optimization, they arc eventually brought to satisfaction within tolerance TOL, indirectly, by virtue of eqs. (14), 
(15), (15a), (15b), and (9b), when the overall procedure converges. 

Modification 2 

An algorithm introduced in reference 10 for the purposes of fitting an empirical function to a set of 
experimental data points may be adapted as a formulation of the subsystem level optimization problem. The 
algorithm requires augmentation of {Y}j by an additional independent variable, LIM X , that also doubles for 


the objective function in a following formulation of the subsystem problem for the ith beam 

“find {Y}j and LIM X such that LIM X is at minimum subject to constraints” (16) 

9k{Y}i <0, fc = l ...NGBi (16a) 

—LI Mi < F Ai < +LIMi (16bl) 

-LIMi < Fa < -VIA M { (16b2) 

{YL}, < {Y}< < {YU}, (16c) 

LIMi > 0 (16d) 
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rod.Tr^tn Tni ^ ? 0l " t,O f n produces the values of UM itnia and { Y} iopt . The LIM iaia is eventually 
reduced to TO I, owing to the following changes in the system level optimization of cq. (9): deleting eq (9c) 
as in Modification 1 above, and replacing C { with LIM { which changes cqs. (9b) and (10) to ' 


LI < TOL, i = 1 . . . NE 


(17) 


uM <m - < uuiu* + % £ ^ {Xj - Xjo) , 

i= 1 • • • NE, j = 1 . • - NX, q = 1 . . . NPi 


(18) 


As in Modification 1, the optimum sensitivity analysis that follows the solution of cqs. (15), (15a) and 
x/S b ^. Carr o icd L 0ut U3in S thc a lg° r *thm from reference 3 instead of reference 2. Neither Modification 1 
previously Cat '° n ^ anythin S m thc step-by-step prescription for thc overall procedure OLD described 

Extension to a Generic System 

The above algorithm lends itself to a complete generalization by replacing the variables and terms specific 
to the framework example with their generic counterparts, while leaving thc organization of the decomposed 
optimization and thc formulations of its elements unchanged. This section explains the substitution of 

the variables and terms in each of the system and subsystem level analyses, subsystem and system level 
optimizations, and thc overall procedure. 

System Level Analysis 

System Analysis, designated SA, is regarded as a black box that converts an input vector {SI} of length 
NSI into an output vector {50} of length NSO. Vector {57} contains as a subset thc vector {AT} of length NX 

that comprises thc system level design variables. Vector {50} also contains {X} as a subset (passed through 
SA). By definition of SA v 6 

50 = /(57) (19) 

assumed differentiable up to thc first derivatives. Consequently, a Jacobian matrix of the first derivatives exists 


050 

051 


dSOj 

~ [dS/j-J ’ l ~ 


1 ... NSO, j = 1...NSI 


( 20 ) 


Equations (19) and (20) form a functional statement that is a very general one and includes two important 
special eases: the first derivatives in thc above Jacobian, cq. (20), degenerate to zero wherever a particular 

50^ ,s not influenced by a particular Sly, and thc derivatives default to unity wherever, for a particular pair 
ij, there is SO( = Slj = Xfc. 

The Jacobian matrix, cq. (20), is obtainable by finite differencing on 5A or by quasi-analytical sensitivity 

ana ysis embedded in SA. In either ease the Jacobian matrix is regarded as an additional output from SA 
separate from SO. 1 


e J hc ' nput dcfincd ln tablc 1 is an example of SI, while thc subset of A’s and 7’s in that input is an example 
of X. The output defined in that tabic is an example of SO, and thc A’s and 7’s present in that output are 
examples of X as a subset of SO that is being passed through SA unchanged. The derivatives dN/dA defined 
in table 1 constitute an example of thc elements of dSO/dSI. 

Subsystem Level Analysis 

It is assumed that there arc NSS black boxes representing thc next lower level of subsystems. As daughters 
o thc SA parent black box, they receive their input in part from thc system level and, in part, from the outside 
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world, but not from each other. Therefore, the subsystem black boxes arc mutually independent and may be 
executed concurrently. The ith black box is designated SSAi. It converts the input vector {SSI}: of length 

A 55/* into an output vector {550}* of length /V550*. 

Vector {55/}, contains selected elements of {50}, some of which may be the elements of X. It includes 
also the vector of the subsystem level design variables {Y}{. 

At this point, it is convenient to define a vector {P}„ * = 1 . . . NP h as a subset of those elements of {55/} • 

which arc functions of X. Naturally, this includes the elements of X, if any arc present in {SSI}:. An example 
of such a subset is cq. (4). J r 

We now assume that for the elements of {P}„ the {55 A}, contains a functional relation, 

{**(*)}< = ( 21 ) 

An example of the above relation is illustrated by cq. (1). However, cq. (20) is more general than eq. (1) 
since it recognizes that not only the elements of X but also other data in 55/*, selected from SO mav be 
computable as functions of {V},. ’ 

For the ith subsystem optimization, it will be useful to define a “discrepancy” vector { DJ5}* defined as 

i DIS jh = ({ Pj)i ~ f({Y)i))/{Pj}i, j = l...NP { (22) 

Examples DISj arc FA, and FI t in cq. (13). The 55A, output {550}, contains {D/5}* and other behavior 
variables of interest. 

Subsystem Level Optimization 

SSOPTi is an optimization of subsystem i. It can be defined by either one of the two alternative ways 
consistent with Modification 1 (cqs. (15), (15a), and (15b)) or Modification 2 (eqs. (16), (16a-A)). 

In the first alternative, corresponding to cqs. (15), (15a), and (15b), the design variables are the elements 
°f the objective function is 

Ci({Y}i, TO = KS({DISj}i,—{DISj}i) 

= {l/p)\n ^£3 ( CX P (pi DI s j}i) + exp (- p{DISj }*)) j , j = 1 . . . NPi (23) 

In this alternative the optimization problem comprises the above C* as the objective function and its 
formulation is as follows: 

find { Y}{ such that {P}i) is at minimum subject to constraints” (24) 

5Jt({^}i) < 0, k = 1 . . . NGSSi (24a) 

{YL}i < {V}* < {YU)i (24b) 

In the second alternative, consistent with cq. (16), the design variables arc the elements of {V}* and an 
additional variable L/M*. The latter doubles for an objective function so that the optimization problem is 

find {V}* and L/A4* such that LI A/* is at minimum subject to constraints” (25) 

< 0, k = 1 . . . NGSSi (25a) 
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-LI Mi < {DiS^i < +LlMi , j = 1 ...NPi 
{YL)i < {Y}i < {YU)i 
LI Mi > 0 


(25b) 

(25c) 

(25d) 


In both cqs. (24) and (25), the constraint functions g k arc evaluated using the data in SSO*. The SSOPTi 
output is {y} iopt an d, dcpcndcntly on the choice of cqs. (24) or (25), F* min = C* m j n or F imin = LIM imin , 
respectively. 

Optimum Sensitivity Analysis and Extrapolation of the Minimum of the Objective Function 

The Optimum Sensitivity Analysis, OSA*, uses the algorithm of reference 3, applied as a post-processor to 
SSOPTi , to yield for each subsystem i a vector of the derivatives of the minimum objective function F* ^ with 
respect to the parameters P. The vector of these derivatives, has the length NP*. The definition of 

P (see discussion preceding cq. (21)) and the functional relationship defined by cqs. (19) and (20) imply that 

{Pj}i = /(X*), j = 1 . . . NPj, k=l...NX (26) 


and that the derivatives exist. These derivatives default to unity for those pairs jk for which Pj = X k . 

Using the above, one may extrapolate Fi by means of the chain rule as an approximate function of X 


v — > r — > 0 1 4 i 

F i = + ~Q~p 

j Q f 


OR 


„ dXj 


P-( Xj-X j0 ), i = 1 . . . NSS, j = l...NX, q = l...NPi 


(27) 


where ( Xj) Q is the X for which the system and subsystem analyses were carried out prior to the subsystem 
optimization. For an example of cq. (27), see cqs. (10) or (18). 

System Level Optimization 

Equating F* either to C* computed from cq. (23), or to LI Mi defined for cq. (25), dependency on the 
choice of Modification 1 (cq. (15)) or Modification 2 (cq. (16)) for the SSOPT formulation, the system level 
optimization, designated SOFT , may be formulated so as to accommodate both alternatives. Defining an 
objective function FS(X) and the vector of constraints G'(X), both computed from the elements of SO, the 
SOPT formulation is as follows: 


“find X such that FS(X) is at minimum subject to constraints” (28) 

G k {X) < 0, k = 1 . . . NGS (28a) 

Fi{X) < TOL, i = 1 . . . NSS (28b) 

XL <X< XU (28c) 

consistent with cq. (9) modified as described in the discussions of the Modifications 1 and 2. 


In the above optimization problem, Fj in cq. (28b) is extrapolated by cq. (27), analogous to the extrapolation 
of Ci in cq. (10) or LIMi in cq. (18). 

Initialization 

It is a recommended practice to initialize the entire optimization procedure by first setting the values of 
{Y}i, i = 1... NSS, and then computing X from cq. (21), recalling that, by definition X is a subset of F. 
This guarantees starting in SA at the system level with the X values that arc physically realizable in SSOPT* 
by {y}* within the YL and YU limits. This operation is abbreviated I NIT. 
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Generic Two-Level Optimization Procedure 

The procedure is the same aa the one described previously for the framework example, restated in the 
generic terms defined in this section. 

1. Execute INIT ; 

2. Execute SA\ 

3. Execute SSOF x f followed by OSA^ i — 1 . . . NSS\ concurrently, if desired and if the computing 
equipment permits; 

4. Execute SOPT ; 

5. If termination criteria set by user arc not satisfied, reset A and Y to A 0 pt and and repeat from 
step 2, otherwise EXIT. 

The procedure output is the optimal data for A, Y , C(A), <;(A), and FS( A). Owing to the formulation 
of SSOPTj in cq. (23) or (24), the optimal values of A and Y will satisfy eq. (20) within TOL. The procedure 
flowchart is depicted in figure 2. 

Concluding Remarks 

Two new techniques arc presented for coupling the levels in optimization by decomposition. The 
techniques constitute improvements of a previously published algorithm for two-level Optimization by Linear 
Decomposition (OLD). The OLD algorithm has been summarized and illustrated by an application example 
to show how the new techniques arc implemented by local modifications in that algorithm. The resulting two 
alternative formulations improve OLD by removing the potential for numerical difficulties that occasionally 
were caused in the original algorithm by an explicit handling of the equality constraints which constituted 
the key coupling between optimizations at two levels. Both alternative formulations eliminate the explicit 
presence of these constraints while satisfying them indirectly. The alternative formulations allow the use of a 
variant of the optimum sensitivity analysis that docs not require second derivatives of behavior and, therefore, 
is computationally less expensive than the variant used in the original algorithm. Thus, an additional benefit 
expected from the modified algorithm is a reduction of its overall computational cost. It is shown that the 
improved algorithm may be generalized to multidisciplinary system applications. 
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Tabic 1. Input and Output of Assembled Framework (System Level) Analysis 


Auh 


{ SSC } 

{Q) 

E 


Input 

cross-sectional area and moment of inertia for ith beam (i = 1 ... N E 1 the 
number of elements in the example is NE = 3), both assumed constant along 
the beam length 

vector of the structure shape coordinates in a reference coordinate system that 
defines locations of the frame support and corner points 

vector of loads Qj , j - 1 ... IS! DOF , applied coincident with the structure 
unsupported degrees of freedom whose number is N DOF 

the material Young’s modulus 


Output 

{{/} vector of displacements t/j., k = 1 . . . N DOF 

{N)i, i - 1 . . . NE vector of the end-forces on the beam, at each end there are three such forces: 
axial force, transverse force, and bending moment 

Lj length of the ith beam 

j4j, 1 ^ and E as an input passed through the analysis to output; 

Derivatives: D(U,A), D(U,I), D(U,SSC), D(N,A), D(N,I), D(N,SSC), D(L,SSC) 


Table 2. Input and Output of ith Beam (Subsystem Level) Analysis 

Input 

NAL and SAL normal and shear allowable stresses, respectively 

Input selected from the Output in table 1: 

vector of the end-forces on the beam 
length of the ith beam 

vector of the beam cross-sectional dimensions Y^, k = 1 . . . NYj, shown in the 
inset, Fig. 1. 

Output 

vectors of the normal and shear stresses, respectively, at judiciously chosen 
points on the end cross sections 

beam displacements in the beam local coordinate system 

{NCR}i and {SCR}i vectors of the normal and shear critical stresses, respectively, for evaluation of 

the beam column and local buckling constraints 


Ni, i = l...NE 

u 

Ai, Ii, and E; 

{Y}i 

{SN}i and {SS}i 

{UB}i 
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Fig.l A portal framework. 






Ln 


Fig.2 Flowchart of the generic procedure. 
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